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1. Introduction 



Advances in both computer capabilities and numerical methods now make it possible to study 
physically light pions in large spatial volumes using a chiral lattice fermion formulation. This 
allows the study of standard quantities such as particle masses and decay constants with enhanced 
precision, removing the uncertainties associated with using chiral perturbation theory to extrapolate 
to physical quark masses. This ability to work directly with physical quark masses also allows us 
to tackle the calculation of more complex quantities where the effects of using unphysically large 
quark masses may be more difficult to estimate. In this paper we will discuss the calculation of 
three such quantities, the decay amplitudes Aq and A2 for the decay of a K meson into 1 = and 
1 = 2 two-pion final states and the mass difference, AM% between the Kl and K$ neutral kaons. 

For K^rTiTi decay we will discuss the calculation of the AI = 3/2 amplitude, A2, with the kaon 
and pion masses and the pion relative momenta taking their physical values [jl|, ||]. By working at 
physical kinematics, the dominant errors come from computational issues such as non-zero lattice 
spacing and finite volume which can be reduced in future calculations by simply working at smaller 
lattice spacing and larger volume. We do not need to deal with an uncertain theoretical framework 
to correct for the absence of the physical % — % relative momentum or unphysical pion or kaon 
masses. For the more difficult calculation of the AI = 1/2 amplitude, Aq and AMk, where such a 
physical calculation is not yet possible, we concentrate on developing the computational methods 
which should allow a calculation with physical kinematics in the not-too-distant future. 



2. K^tnn decay 

Since the weak interaction W ± bosons which mediate this decay are far too massive to simulate 
in a lattice QCD calculation, our first step must be to represent W ± exchange by the effective four- 
Fermi interaction which results if the W exchange process is treated as taking place at a space-time 
point. The resulting effective Hamiltonian is written as 



= 1) 




/ \ V t d v t * 



a > . (2.1) 



Here V qq > is the Cabibbo-Kobayashi-Maskawa matrix element connecting the charge —1/3 quark 
q' to the charge +2/3 quark q. The Wilson coefficient functions y,-(ju) and Z/(jU) depend on the 
scale /i at which the four-quark operators <2, are normalized and have been determined in QCD and 
electro-weak perturbation theory through second order; see Ref. [^] for a thorough discussion. 

The ten dimension-six, four-quark operators Qj, 1 < i < 10 are not independent but arise 
from particular phenomena and are defined in Eqs. (4-23) of Ref. The operators Qi, i = 1,2 
represent the naive current-current interaction resulting from simple W exchange and transform 
under SUl(3) x SUr(3) as both (8,1) and (27,1) representations. Those with i = 3 through 6 
result from QCD penguin graphs in which a quark-anti-quark pair emerges from the point-like W 
exchange and annihilate into a gluon which subsequently creates a possibly different q — q pair. 
These operators transform in the (8,1) representation. The final four operators arise from electro - 
weak penguin diagrams in which a quark-anti-quark pair emerging from the point-like W exchange 
annihilate to create a photon or Z boson and transform as (8, 8), (8, 1) and (27, 1). While one order 
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smaller in electro-weak perturbation theory, the (8,8) operators can be important because of the 
suppression of the leading order, AI = 3/2 component of the usual (27, 1) matrix elements which 
arises from the "A/ = 1/2 rule" and an extra power of M\ appearing in chiral perturbation theory. 

While this historical classification of the ten <2; operators is convenient for discussing the 
electro-weak phenomena underlying the decay, when discussing renormalization and operator mix- 
ing it is better to work with a set of seven linearly independent operators, four transforming in the 
(8, 1) representation, one in the (27, 1) and two transforming in the (8,8) representation. If renor- 
malization conditions are imposed in which the flavor-symmetry breaking quark masses can be 
neglected then these three classes of operators will renormalize independently. 

Critical to an accurate lattice calculation of a weak decay process such as K — > %% is the ability 
to relate the lattice-regulated operators appearing in such a calculation to the continuum operators 
for which the Wilson coefficients were originally computed. This can be done with increasing pre- 
cision using the intermediate Rome-Southampton RI/MOM approach [j5|, enhanced by a number of 
refinements over the past decade. Here one introduces a regularization-independent scheme to nor- 
malize the lattice operators in which particular Landau-gauge-fixed, spin-color-projected, Green's 
functions are normalized at large, off-shell momenta at a scale characterized by /I. In this way, the 
use of lattice perturbation theory is avoided and the original bare lattice operators are expressed in 
terms of operators renormalized in a scheme which has a well-defined continuum limit. Because 
of this non-perturbative step, these methods are referred to as non-perturbative renormalization 
(NPR). If jU is sufficiently large, these same conditions can be accurately imposed in a perturbative 
calculation allowing this RI/MOM renormalization to be connected to the standard perturbative 
MS continuum scheme in which the Wilson coefficients are typically evaluated. 

The effectiveness of these NPR techniques can be seen in the size of the normalization er- 
rors presented below. Three important developments, all used in the calculations discussed here, 
have made these NPR techniques more accurate. The first is the use of non-exceptional momenta 
when imposing the RI/MOM normalization conditions [g]. This choice of momenta makes these 
normalization conditions infrared safe, significantly reducing the contributions of small loop mo- 
menta and making the use of QCD perturbation theory to relate the RI/MOM and MS schemes 
more reliable. The second improvement recognizes that even at 0(a 2 ), lattice artifacts appear in 
the RI/MOM normalization conditions which are not 0(4) invariant. As a result the amplitudes 
being studied depend not only on the usual Lorentz scalars such as p 2 or p ■ p', where p and p' are 
lattice momenta, but also on the direction of these momenta relative to the 4-D lattice axes. This 
few-percent direction dependence introduces irregularities into what should be smooth dependence 
on p 2 and prevents the evaluation of the continuum limit if an a 2 extrapolation is attempted using 
lattice momenta in different directions and consequently with different a 2 corrections. These diffi- 
culties can be avoided by imposing twisted boundaiy conditions on the fermion propagators used 
to evaluate the RI/MOM normalization conditions. By vaiying the degree of twist (the fermion 
phase change when passing through a boundary) the magnitude of the lattice momentum can be 
varied, e.g. to compensate for a change in lattice spacing when extrapolating to the continuum 
limit, without changing the direction of that momentum relative to the underlying lattice [|7| [8j|. 

The third advance in NPR methods is using the momentum dependence of the RI/MOM renor- 
malization factors to relate operators renormalized at significantly different scales [J7J [8]]. This al- 
lows operators which are evaluated in a coarse lattice calculation, such as the K — > %% amplitudes 
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discussed below to be accurately renormalized. In this K — > %% calculation the physical size of 
the lattice is increased by using a relatively large lattice spacing of 0.144 fm. This limits the size 
of the external momenta that can be employed in an RI/MOM renormalization condition to < 1 
GeV, a scale too low for the application of perturbation theory. However, by using the step-scaling 
function determined in a companion calculation on a finer lattice, these operators renormalized at 1 
GeV can be accurately renormalized also at 3 GeV, a scale at which a perturbative calculation can 
be used to convert to the standard MS scheme. 

The final difficulty that must be overcome to carry out a realistic lattice calculation of K — > Tin 
decay arises from the non-zero relative momentum of the physical final state pions. In the usual 
lattice QCD calculation, one uses the Euclidean time-development operator e~ Ht at large t to project 
onto the QCD energy eigenstates with the lowest energies. While this strategy works well to 
construct the initial K meson, when applied to a state of two pions, unphysical, threshold states 
with zero relative momentum pions result. 

This difficulty can be overcome in two steps. The first, proposed by Lellouch and Luscher [Q], 
recognizes that in finite volume there will be additional excited states containing two pions obeying 
energy quantization conditions that depend on the volume. With an appropriately chosen volume, 
the resulting two-pion energy can be adjusted to equal that of the decaying kaon. Thus, by identi- 
fying the transition amplitude to this excited state and introducing a finite- volume correction factor 
derived in Ref. ^ the physical, on-shell, K — > %% decay amplitude can be computed with con- 
trolled errors. However, the extraction of such an excited state is typically difficult and a second 
step of imposing boundary conditions to remove some or all of the energy-non-conserving, lower 
energy % — % states from the calculation can substantially improve the result [10, 11, 12]. We will 
discuss such boundary conditions in greater detail below. 



2.1 K^titi decay with AI = 3/2 

The K — >• 7171 decay amplitude most accessible to the methods of lattice QCD is A2 which 
describes the decay into the n — 71 state with isospin 2. This state does not have vacuum quantum 
numbers and quark flavor conservation implies that the valence quark lines connect the initial kaon, 
the effective four-Fermi weak operator and the final two pions. Thus, there are no disconnected 
diagrams, no need for a vacuum subtraction and the light sea quarks do not play a critical role, 
allowing boundary conditions to be applied only to the valence quarks with resulting errors that 
vanish exponentially in the lattice size JTH]. 

We ensure that the / = 2, % — % final state in a physical K decay is also the lowest energy 1 = 2 
state in our lattice calculation by imposing anti-periodic boundaiy conditions on one of the quarks 
in each pion and adjusting the lattice volume so that n/L is close to the 205 MeV momentum of 
the physical decay pions. This is accomplished in two steps. First we use isospin symmetry to 
relate a physical AI = 3/2, K decay amplitude to a transition caused by a related weak operator 
carrying charge +1: K + — > 7l + n + . We then impose anti -periodic boundaiy conditions on the d 
quark. The resulting K + meson will obey periodic boundary condition so the initial K + can be 
given zero momentum. However, each 7i + mesons will have non-zero momentum equal to 71 /L 
except for rescattering effects. Normally the use of such isospin violating boundary conditions 
would introduce potentially dangerous mixing between the 1 = 2 final state of interest and the 1 = 
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state which overlaps with the vacuum. However, by arranging the final state to carry charge 2 we 
have assured that no / = final state is possible. 



The first calculation adopting this strategy was C. Kim's Ph.D. thesis Q14[ ] which used a 
quenched 16 3 x 32 ensemble and unphysically massive kaons and pions. An interesting related 



calculation [ ]15[ ] studied the decay of a kaon carrying non-zero momentum, an alternative technique 
to achieve a final % — % final state with physical relative momenta. The statistical errors and re- 
sulting practical restriction on the separation between the kaon and % — % sources in this non-zero 
momentum approach appear to favor the boundary condition strategy discussed here. 

We now turn to the recent RBC/UKQCD calculation of A2, using a large lattice spacing of 
l/a = 1.364(9) GeV, a near physical (partially quenched) pion mass of 142.11(94) MeV, a kaon 
mass of 505.5(3.4) MeV and a two-pion energy of 484.5 (4.2) MeV [jl|, 0]. This two-pion energy 
corresponds to a d quark which obeys anti-periodic boundary conditions in two of the possible three 
spatial directions. Three independent effective weak operators contribute to the complex amplitude 
A2' Q( 27,1 \ Q( s -^ and gt 8 8 )'". The computational setup is shown schematically Fig. 1(a) . The 
light and strange quark propagators are computed using both periodic and anti-periodic in the time 
direction. The strange quark source is located at the tg, where seven values of tx are used varying 
between t % + 20 and t n + 44 in steps of four time units where t % locates of the source of the light 
quarks. The sum of the periodic and anti-periodic propagators are used for both the light and 
strange quarks to increase statistics and reduce around-the-world effects. 




O This Calculation 
□ Hoogland et. al. (A) 
o Hoogland et. al. (B) 
a Lostyet.al. 
— ■ Phenom. (Schenk) 



k (MeV) 



(a) The setup for the calculation of A2 with the (b) Comparison of calculated phase shift 
kaon wall source at the time the tg, the weak op- with experimental results [F^, |n|] and a 
erator at t and the two-pion source the time t n . phenomenological curve Q18|]. 

Figure 1: 

In Fig. ^ we show a ratio R l {tQ) of three-point to two-point correlation functions from which 
the lattice matrix elements of the three operators g( 271 ) , g( 8 > 8 ) and Q^^ m are determined: 



(2.2) 



Here the location of the two-pion source has been set to t % = and the label i distinguishes the three 
Al = 3/2 operators being studied. The amplitudes ^ are the lattice matrix elements we are trying 
to determine while Zg and Z %% are source normalization factors which can be determined directly 
from the two-point functions. The ratio /?'(?g) should not depend on tQ if only kaon and two-pion 
states are present. The solid line shown in each graph indicates the fitted results, the dotted lines 
the width of the error band and the horizontal position of these lines, the range over which the fit is 
performed. 
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Figure 2: Ratios of three point and two point functions described in the text from which the matrix elements 
of the three operators g( 271 ) , g( 8 - 8 ) and g' 8 ' 8 '™ are determined, plotted as a function of the time tg at which 
the operator is located. These are computed using anti -periodic boundary conditions in two spatial directions 
for our full ensemble of 146 configurations. The two-pion source is located at t K = and the kaon source at 
t K = 24. 



The complex amplitude A 2 can be obtained from the lattice matrix elements ^ using: 

yfi^K,lQ(M)zo-4 /2 



2%q n U dq % dq n 



(2.3) 



obtained from a combination of Eqs. (18) and (20) in Ref. [0J. The Wilison coefficients Q(/i), 
evaluated at \i = 3 GeV are obtained from next-leading-order formulae in Ref. [Qj. The renor- 
malization matrix, Zy, transforms the lattice operators used on the DSDR ensemble into the MS 
scheme at \i = 3 GeV. This matrix is determined using the step scaling methods described above. 

The square brackets contain the Lellouch-Luscher correction factor with its sum of derivatives 
of a known kinematic function and the I = 2 n — n s- wave phase shift 8. This phase shift can 
be determined directly by applying the Luscher quantization condition [19] to the calculated n — % 
energies for zero and two twists. The results are shown as the two red open circles in Fig. 1A and 
contribute only 6% to the factor in square brackets. 

Table [I] lists the systematic errors in this calculation. The largest comes from the large finite 
lattice spacing and is estimated in two ways. First we vary the physical quantity that is used to 



determine the lattice scale, using ma, fn, fie and tq. Since the lattice spacing enters Eq. (2.3) with 
the third power, this explicit, «5% uncertainty is amplified 3x. We obtain a similar estimate from 
the a 2 dependence of the K°—K° matrix element of the related (27,1) operator which determines Bk 
and has been computed at a number of lattice spacings. A second error that should be mentioned 
arises from partial quenching of the light quark mass, the valence and dynamical pions having 
masses of 142 MeV and 171 MeV respectively. We estimate an upper bound on the resulting error 
from threshold calculations of A2 performed on the RBC/UKQCD l/a = 2.28 GeV ensembles in 
which the dependence on the sea quark mass was zero within statistical errors. 

We find Re(A 2 ) = 1.381(46) stat (258) syst lO" 8 GeVandIm(A 2 ) = -6.54(46) stat (120) syst l(r 13 
GeV Here ReA 2 agrees well with the experimental values of 1.479(4) and 1.573(57) 1(T 8 GeV 
obtained from K + and K$ decays respectively. The difference between these two experimental 
numbers results from isospin breaking effects, which are not included in our calculation. The 
imaginary part of A 2 is unknown so that this result represents its first direct determination. 
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ReA 2 


ImA 2 


lattice artifacts 


15% 


15% 


finite-volume corrections 


6.0% 


6.5% 


partial quenching 


3.5% 


1.7% 


renormalization 


1.8% 


5.6% 


unphysical kinematics 


0.4% 


0.8% 


derivative of the phase shift 


0.97% 


0.97% 


Wilson coefficients 


6.6% 


6.6% 


Total 


18% 


19% 



Table 1: Estimates of the major systematic errors in this calculation of Re A 2 and YmAz- 




type:! 



ty|->e4 



Figure 3: The topologies distinguishing the four types of diagram contributing to the 7 = amplitude Aq 



2.2 K-^kk decay with AI = 1/2 

The calculation of the AI = 1/2 amplitude Aq describing kaon decay into the I = 0, % — % state 
is much more difficult than that for A 2 . A total of 50 different contractions contribute which can 
be organized into the four types shown in Fig. |3| The greatest difficulty is cause by disconnected 
diagrams shown as type 4. Such diagrams lead to a signal-to-noise ratio which decreases expo- 
nentially with increasing time separation and imply that both sea and valence quarks enter in the 
physical propagating states. This requires that if boundary conditions are used to remove unwanted 
zero relative momentum % — % states, these conditions must be imposed both when computing the 
valance propagators and when generating the gauge ensembles. A further difficult comes from the 
quadratically divergent quark loops found in diagrams of type 3 and type 4. While these terms do 
not contribute to on-shell matrix elements, they can enhance off-shell, excited state contributions 
by factors of 10-20 and some partial subtraction must be carried out if the usual large-time methods 
are to be able to successfully remove the resulting excited state contamination. 

However, while severe, these difficulties may be more easily overcome for a quantity such as 
Aq which involves light pions with their positive definite propagators and a kaon, which has half the 
mass of the much more difficult nucleon. Here we will summarize recent results for Aq obtained 



with unphysical, threshold kinematics and relatively heavy pions which appear in Ref. [gOJ and 
in the Ph.D. thesis of Q. Liu [pj]. The former were obtained from 800 configurations using 2+1 
flavors, l/a = 1.73 GeV and an 16 3 x 32 lattice volume. In contrast to the calculations of A 2 
described earlier in which the two-pion was fixed on a single time slice, the correlation functions 
used in the calculation of Aq were computed for each of the possible 32 time slices. With these large 
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(a) Q2 , which gives the dominant contribution to (b) 26 » which the dominant contribution to Im(Ao) 
Re(A ) 

Figure 4: The point functions for the operators Q2 and as a functino of t the time slice over which that 
operator is summed. In this figure, the pions in the two-pion source are located at the times —4 and and 
the kaon source at time 12. The close points show the result if the disconnected diagrams are omitted while 
the open point show the complete result. 



statistics, it was possible for the first time to determine Re(Ao) from an explicit K — > %% matrix 
element, with a statistical error of « 25%. Of course, this 25% error was achieved for unphysical 
threshold kinematics in which the two pions are at rest and for an unphysically large pion mass 
of 422 MeV. (The valence strange quark mass was adjusted to make m# = E nn so that an energy 
conserving decay was studied.) The imaginary part of Ao could not be distinguished from noise. 
This quantity is dominated by QCD penguin diagrams of type 3 and is more difficult to compute. 

Exploiting what was learned in this 16 3 x 32 calculation, a more ambitious 24 3 x 64 calculation 
was undertaken [21], again using the same gauge action and a smaller light quark mass, giving a 329 
MeV pion mass. Two important improvements were realized. First, propagators both periodic and 
anti-periodic in the time were evaluated. By using the sum of these propagators, the distance to the 
nearest periodic image of the source is moved from 64 to 128 time units, substantially reducing the 
around-the-world effects which gave the largest excited state contamination from the divergent type 
3 diagrams in the 16 3 x 32 calculation. Second, the two-pion source was modified so that the pions 
were emitted from different time slices. (We refer to this as a split-pion source.) This substantially 
reduced the coupling to the vacuum state and the resulting noise in the disconnected diagrams. 
With these advances, improved results were obtain using only 137 instead of 800 configurations. 
Now both the real and imaginaty parts of Ao can be resolved. The three-point functions for the 
two dominant matrix elements are shown in Fig. [| The open symbols show results without the 
disconnected graphs while the closed symbols show the full result. While the statistical errors are 
much larger for the full amplitudes, the central values are consistent with those coming from only 
connected graphs. At present the disconnected diagrams are only a source of noise. 

With the demonstration of methods capable of resolving both the real and imaginary parts 
of Ao, we can now work toward the physical kinematics used for the calculation of A2 described 
above. The same Iwasaki plus DSDR gauge action but with the sea quark mass reduced to its phys- 
ical value should permit calculation with physical pion and kaon masses. Recent RBC/UKQCD 
experience [ 22 ] with the use of the Mobius variant of DWF [ 23 , ^] suggests that the extent in the 
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fifth dimension of 32 used in the current Iwasaki + DSDR, DWF calculations can be reduced to 16 
to enable a high statistics study of Aq. More challenging is the needed non-zero relative momentum 
of the two final-state pions. As discussed above we plan to accomplish this by imposing G-parity 
boundary conditions on the sea and valence quarks. This effort is well underway in a project carried 
out by C. Kelly. Quenched studies have been performed [25] and code written which is currently 
under test to generate Nf ■ = 2 + 1 gauge configurations obeying G-parity boundary conditions. 

While a large reduction in the errors from the disconnected diagrams results from the use of 
split-pion, Coulomb gauge fixed wall sources, we are optimistic that all-to-all propagator meth- 
ods [ p6| , 27 ] will allow us to construct even more effective split-pion sources using localized pion 
wave functions, possibly gaining a further factor of two reduction in statistical error. A final criti- 
cal acceleration is provided by the advance of computer technology. With the now available BG/Q 
computer hardware and highly efficient QCD code described in Peter Boyle's talk at this meet- 
ing [28], the generation of 3K time units of a specialized 32 3 x 64 G-parity ensemble requires only 
a few months on a BG/Q 1024-node rack. The difficult, high-statistics measurements used to obtain 
the 24 3 x 64 results for Aq presented here formed the final project carried out on a 4096-node QC- 
DOC partition in 2011. The sustained performance and memory size of this partition (1 Tfl ops/500 
Gbytes) are equaled by 32 BG/Q nodes, allowing a calculation on a 512 node BG/Q partition to run 
16 times faster. While important aspects of a calculation of Ao (which will give the standard model 
prediction for e' /e) with physical kinematics will remain uncertain until large-scale experiments 
are begun in a few months, this calculation should be possible with present resources. 



3. Computing the K L — K$ mass difference 

Much must yet be accomplished to accurately carry out the K — > %K decay calculations dis- 
cussed in the previous section. However, it is plausible that the basic methods are now understood 
and that presently available measurement algorithms and computer resources will be sufficient for 
the task. In this last section of the talk, we will discuss a more difficult topic at a much earlier stage 
in development, the calculation of "long distance" contributions to second order weak processes. 
One might expect a large difference in complexity between the first order (one-W exchange) pro- 
cesses responsible for K — > Tin decay and the second order weak (two-W exchange) processes 
needed for K° -K° mixing. However, if the two-W exchange is dominated by momenta on the 
order of the W boson or top quark mass, then at the scale of hadronic phenomena this second order 
process will appear to take place at a point and can be represented by a four-quark operator in a 
lattice QCD calculation just as is done for processes that involve a single W exchange. 

The CP violating, K — K° mixing amplitude is dominated by short distance effects and the 
corresponding mixing parameter Ek is typically expressed as a product of a short distance Wilson 
coefficient, derived from a box diagram in which two W bosons are exchanged, and the low energy 
matrix element of a four-quark operator evaluated between K° and K° states, which determines the 
familiar Bk parameter. However, there is also an « 5% "long distance" contribution to Sk [ |29| ] 
in which the two exchanged W mesons are separated by distances on the order of 1 /Aqcd or 
\jm n . Further, the CP conserving part of the K — K° mixing amplitude which gives the Kl — K$ 
mass difference AM% , receives a potentially large, long distance contribution. One expects that the 
largest contribution to AM K comes from distances on the order of the inverse charm quark mass, 
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Figure 5: The two possible low energy descriptions of a box diagram entering K° — K° mixing. Note there 
can be overlap between these long- and short-distance contributions. The long-distance part, expressed as 
the product of two weak operators Hw, will itself contain a short distance part when these two operators 
collide in space-time. Such a short distance component of this long distance amplitude will typically be 
incorrect but, if significant, can be removed by explicit subtraction. 



l/m c . This is conventionally referred to as a short distance contribution. However, perturbation 
theory calculations show large NNLO terms [|^] suggesting that this scale might be better thought 
of as also a relatively long distance at which perturbation theory has become unreliable and again 
a lattice QCD calculation is needed. These short- and long-distance contributions to a box diagram 
which enters K° — K° mixing are illustrated in Fig. |[ 

Two important questions arise when considering a lattice QCD calculation of long distance 
effects in second-order weak amplitudes. The first is associated with the use of Euclidean space 
methods to compute amplitudes with possible real intermediate states. These states result in prin- 
cipal parts appearing in formula for the real parts and give imaginary parts which cannot appear in 
a Euclidean space calculation. The second issue is the need to develop effective lattice methods to 
evaluate the resulting four-point functions in which a K° is transformed into a K° by the action of 
two, separated, first order weak operators H\y. We will now consider each of these issues in turn. 

The computation of long distance parts of second order weak processes is key to many impor- 
tant weak processes, for example rare K decays @]. Here we will focus here on K° — K° mixing, 



typically described in the Wigner-Weisskopf formalism by the evolution equation [ 32 ] : 




(3.1) 



where the 2x2 matrices M and T are given by: 



27r£ / dE{i\H w \a(E)){a(E)\H w \j)8(E - m K ) 

a J2m n 



M, 



5> 



dE 



2m n 



{i\H w \a(E)){a(E)\H w \j) 
niK — E 



(3.2) 
(3-3) 



We are using the subscripts and to represent the K° and K° states and the generalized sum over 
a and integral over the energy represents the sum over a complete set of energy eigenstates. 
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In the method of Lellouch and Luscher, the Luscher finite volume condition []19|] is used to 



relate the finite volume energy of the K—%% system to the infinite volume % — % scattering phase 
shift evaluated at the resonant energy corresponding to on- shell K — >• %% decay. Both quantities 
are evaluated to first order in H w : the first order shift in the finite volume % — % energy coming 
from the coupling to the degenerate K state and the Breit-Wigner resonant contribution to % — % 
scattering evaluated at the kaon pole. 

A generalization of this approach connects the infinite-volume, second-order mass difference 
AMk which shifts the location of the kaon pole in resonant % — % scattering to the second order 
degenerate perturbation theory of a finite volume system including two or three nearly degenerate 
states: a two-pion state and either a CP even or both CP even and odd K° — ^° states, depending 



on whether we wish to treat only AMk or both AM% and Ek [33, 34]. In contrast to the Lellouch- 
Luscher correction, the effect of finite volume is not a multiplicative factor but an additive, 1/L 3 , 
correction. For AMk the infinite- and finite- volume expressions are related through 0(l/L ) by 



AMr 



21 



(K°\H w \n)(n\H w \K° 



m K - E„ 



+ 



dE 



i a 2 (0+g o ) 2 

2 — dip — \\M H m K s)\ 



dE„ I dE 



(n\H w \K. 



E=E, 



"o 



(3.4) 



Here the state \n ) is a two-pion state whose energy has been adjusted to equal that of the K meson 
by the choice of volume. This result suggests that just as for the K — > %% decay amplitude, the long 
distance parts of second-order weak ^° — K° mixing amplitude should be accessible to a Euclidean 
space calculation with controlled finite-volume errors. 

We now discuss the methods needed to cany out such a second-order weak calculation using 
lattice QCD. This has been an active project of the RBC/UKQCD collaboration for the past two 
years and has been carried out and reported on at both the current and previous Lattice Field Theory 
Symposia by J. Yu j35| , 36]. The basic idea is to integrate the product of two effective weak 
Hamiltonia Hy/(tt), i = 1,2 over a fixed temporal region t a <h,t2<tb, taking the matrix element 
between K° and K° states: 



dt 2 / dhH w {t2)Hw(h)K?{ti) |0). 



(3.5) 



This product is illustrated in Fig. 



We can evaluate Eq. (3.5) by inserting sums over complete sets of intermediate states and 
carrying out the integrals over t\ and t%. These are actually sums over discrete times which can be 
evaluated as partial geometric series and approximated for small a to give: 



Ufco M K~En \ Kh a) M K -E n 

e (M K -E„)(t b -t a )\ , ) 

+ Mk _ e \+-(K°\H w \n )(n \H w \K )(t b -t a ) 2 V 



(3.6) 



This equation contains four terms which should be interpreted: three within the large curved brack- 
ets and the fourth term with the 1/2 prefactor. The first term, proportional to % — t a , is the desired 
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Figure 6: A diagram showing the elements of a lattice calculation of the long distance contribution to 
AMk- The times t\ and to appearing in the two effective weak operators are integrated over a fixed interval 
to < hJi < h- in Minkowski space such an amplitude would yield the second order mass shift times the 
elapsed time interval % — t a . 

cl _ s 

^ — * — \ , ^ — ' \J c, u \ 



(a) Type 1 



(b) Type 2 





(c) Type 3 (d) Type 4 

Figure 7: Examples of the four types of diagram that enter the calculation of the second order weak 
contribution to AMk- 

finite-volume expression of the mass difference. The second is an uninteresting, time-independent 
constant, the third disappears at large time for all states more massive that the kaon. States which 
are lighter than the kaon will give exponentially growing contributions which must be evaluated 
separately and subtracted. The fourth term, proportional to (tf, — t a ) 2 , results if the volume has 
been adjusted to create a state, |«o). degenerate with the K. This term must also be identified and 
discarded if the prescription to control finite volume errors described above is followed. 

Figure shows the four types of diagram which contribute to the amplitude si '. In the first 



exploratory calculation of this quantity [ ]35| , p6| |, diagrams of type 3 and 4 which are disconnected 
in the ( or i channel are neglected. This calculation is performed on an ensemble of 16 3 x 32 
configurations with l/a = 1.73 GeV, generated using 2+1 flavors and the Iwasaki gauge action. 
In order to realize GIM cancellation, a valence charm quark is included and the effective weak 
Hamiltonian appropriate for four flavors is used. This Hamiltonian includes the six weak operators 

Q? =Tif(l- fHqjfil - fWj Qf' = Sit {I ~ 7 5 )^y J (l " T 5 )*?; (3.7) 

where i and j are color indices while q and c( are u and/or c quarks. The introduction of a charm 
quark into a lattice calculation with an inverse lattice spacing l/a = 1.73 GeV will introduce po- 
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(a) The integrated correlator si as a function of 
the integration interval tj, — t a after GIM cancel- 
lation with a 0.954 GeV valence charm quark. 
The red squares and blue diamonds are the re- 
sults before and after the subtraction of the ex- 
ponentially increasing 7t° term respectively. We 
include only the Q\ Q\ operator combination in 
this plot. 



Valence Charm Quark Mass / Gev 

(b) The mass difference AM^) , obtained from 
the slope of the amplitude in Eq. (3.6) with re- 
spect to tj, — 1„ after GIM cancellation and sub- 
traction of the light pion contribution, plotted as 
a function of the valence charm quark mass. The 
'11' superscipt indicates that only the contribu- 
tion of the product of operators of the type Q\ is 
shown. 



Figure 8: Results from the lattice QCD calculation of AMk- 



tentially large discretization errors requiring future work with smaller a. The GIM cancellation is 
complete: the short distance part of our Hw x Hw product is inaccurate at the level of (m c /mw) 2 , 
much smaller than the (m c a) 2 discretization errors. 

Preliminary results are shown in Fig. [| The left panel shows the linear behavior as the integra- 
tion time interval tj, — t a is varied suggesting that the required slope is not difficult to extract, after 
the exponentially growing contribution from the light 71° state has been removed. The right panel 
shows the increasing values for AMk that result for increasing charm quark mass, m c . This figure 
also suggests the presence of a sizable constant term needed to describe the large m c dependence, 
reflecting a significant, m c -independent long-distance contribution, at least for the large values of 
light and strange quark masses used here. The final results for AMk vary between 5.12(24) and 
9.31(65) 10~ 12 MeV as the kaon mass varies between 563 and 834 MeV, for a pion mass of 421 
MeV. These results are somewhat larger than the experimental value of 3.483(6) 10~ 12 MeV. 



4. Conclusion 

Substantial advances in computer capability and powerful new numerical methods are dramat- 
ically increasing the accuracy with which standard quantities can be computed using lattice QCD 
and expanding the range of important quantities which can be calculated. By working at physical 
light quark mass and relatively coarse lattice spacing, it is now possible to directly calculate the 
/ = 2 K — > %% decay amplitude A2 with the largest error coming from finite lattice spacing effects. 
Over the next couple of years, it should be possible to repeat these calculations on a series of ensem- 
bles with varying lattice spacing, reducing the error on A%, perhaps to the 5% level expected from 
isospin breaking effects. The more difficult 1 = amplitude Ao has been computed for unphysical 
kinematics and calculations with physical kinematics are now being actively planned. Much more 
ambitious is the calculation of the Kl — K$ mass difference AMk which appears to be within reach. 
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While exploratory calculations are now underway, results with controlled errors on the 5% level for 
AMk are likely five years away. The results reported here represent a major research direction of 
the RBC and UKQCD collaborations and I thank my collaborators whose work is being described. 
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